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^ : ABSTRACT 

^ I In the light of the recent results of the stellar interferometry, we examine the 

CS| I nature of the extra molecular layer outside the photosphere of red supergiant 

^ I stars, so far studied mostly with the use of the infrared spectra. Although the 

I visibility data are more direct probes of the spatial structure of the outer atmo- 

^ I sphere, it is essential that they are analyzed in combination with the spectral data 

I of a wide spectral coverage. In the case of the M2 supergiant /i Cephei, several 

O ' sets of data, both spectra and visibilities, strongly suggested the presence of an 

\o : 

O ■ 

■ 

Tex ~ 1600 K, column densities of CO and II2O molecules Acoi ~ 3.0 x 10^° cm ^, 
O ' and located at about one stellar radii above the photosphere or ^ 2.0 i?*. 

^ • The result shows reasonable agreement with the one based on the infrared spec- 

■ tra alone, and this may be because the infrared spectra already include some 

■ information on the spatial structure of the outer atmosphere. It is important, 
^ ' however, that the model inferred from the spectra is now fully supported with 

■ the recent visibility data. In the case of the M2 supergiant a Orionis, the infrared 
spectra and visibilities show a consistent picture in that its molsphere is closer 
to the photosphere (i?in ~ 1.3 i?*) with higher gas temperature (Tgx ~ 2250 K) 
and lower gas column density (Acoi ~ 10^°cm~^), compared with that of /i Cep. 
Some controversy on the interpretation of the mid infrared data of a Orionis can 
be reconciled. Given that the presence of the extra molecular layer is reasonably 
well established, consistently with the spectral and visibility data, in at least two 
representative red supergiant stars a Orionis and fi Cephei, the major unsolved 
problem is how to understand the origin of such a rather warm and dense layer 
in the outer atmosphere. 



extra-molecular layer (which we referred to as molsphere for simplicity) , and the 
basic parameters of the molsphere are estimated to be: excitation temperature 



Subject headings: molecular processes - stars: individual (a Ori, /i Cep) - stars: 
supergiants - stars: late-type 
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1. Introduction 

Detailed studies on the structure of the outer atmospheres of cool luminous stars started 
with high resolution infrared spectroscopy. Especially, Fourier Transform Spectroscopy 
(FTS) pioneered by P. & J. Connes (e.g. Connes 1970) and developed at KPNO (Hall et 
al. 1979) played a major role for this purpose (e.g. Ridgway and Brault 1984). The possi- 
ble presence of an extra-molecular layer distinct from the cool expanding wind as well as 
from the hotter chromosphere was first noticed in Mira variables whose spectra of the CO 
2-0 band showed extra absorption that remained stationary against the velocity shift of the 
photospheric lines due to pulsation (Hinkle et al. 1982). In the case of red supergiant stars, 
subtle excess absorption seen on the high resolution FTS spectra of the CO first overtone 
bands was interpreted as due to a quasi-static molecular layer above the photosphere (Tsuji 
1987). The presence of the strong photospheric absorption lines of CO, however, made it 
rather difficult to be fully convincing of such a result. 

By hindsight, more clear evidence for the presence of the extra molecular layers was 
observed with the balloon-borne IR telescope, known as Stratoscope II, at the infancy of the 
infrared astronomy in the early 1960's by Woolf et al. (1964) and by Danielson et al. (1965), 
who have correctly identified water vapor spectra in several early M giant and supergiant 
stars. Later, it was recognized that the Stratoscope data should be more convincing evidence 
for the extra molecular layers, since water cannot be formed in the photospheres of these early 
M giant and supergiant stars (Tsuji 2000a). Also, another space IR observation by Russell 
et al. (1975), which showed fiux excess in the 5-8 /xm region of /x Cephei, was interpreted as 
due to the thermal emission of water in the circumstellar envelope (Tsuji 1978a). The next 
major progress was provided by the Infrared Space Observatory (ISO) mission launched on 
1995 (Kessler et al. 1996), and ISO uncovered that water exists everywhere in the Universe, 
including the early M giant stars (e.g. Tsuji et al. 1997; Decin et al. 2003). Finally, water was 
found in emission in the 6-7 /im region as well as in the 40 /im region of /i Cephei observed 
with Short Wavelength Spectrometer (SWS) (de Graauw et al. 1996) on board the ISO, and 
the extra- photospheric origin of water lines appeared to be more likely (Tsuji 2000b). 

The problem of water in cool stars, however, may involve some other aspects. Especially, 
water vapor lines were detected on the high resolution 12 fj,m spectrum of the K giant star 
a Boo by Ryde et al. (2002), who suggested that this result can be explained as due to an 
anomalous structure of the photosphere. Since it is not hkely that such an early K giant star 
could have an extra molecular layer, it may be difficult to apply the molsphere scenario to 
this case. More recently, water vapor lines were clearly observed in the 12 ^m region of a Ori 
by (Ryde et al. 2006), who argued the possibility of the photospheric origin as in the case 
of a Boo. This observation confirmed the previous detection of the pure rotation lines of 
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H2O in a Ori and a Sco by Jennings & Sada (1998), who suggested that the water hnes may 
be formed in the temperature minimum region of the photosphere-chromosphere transition. 
Thus, to understand the full meaning of water spectra in cool stars (but not necessarily 
very cool such as Miras and IR stars), more works on cool stellar atmosphere, including the 
photosphere, chromosphere, and outer molecular layer, should be required. 

For more detailed analysis of the structure of the extra molecular layers, direct infor- 
mation on the spatial structure should be indispensable. Direct evidence for the presence of 
the extended molecular envelope was first given for Mira variables o Ceti and R Leo with the 
Speckle intcrfcromctry by Labcyrie et al. (1977), who showed that the angular diameters in 
the region strongly blanketed by the TiO bands are larger by two times than those in the 
region free from the TiO bands. In normal M giants and especially in M supergiants, the 
angular diameters in the TiO bands relative to those outside of the TiO bands were found 
to be larger by as much as 20 % by Quirrenbach et al. (1993), who pointed out that such 
an effect should not necessarily be related to the stellar pulsation and hence may be a more 
fundamental property of cool luminous stars. Recently, multi-wavelength measurements of 
the angular diameters are extended to the infrared region, including the K band (Dyck et 
al. 1992, 1996, 1998; Perrin et al. 2004a), V band (Chagnon et al. 2002; Mennesson et al. 
2002), and mid- infrared (Weiner et al. 2003). For the case of a Ori, the resulting apparent 
diameter at 11 /xm is ^ 30% larger than that at the K band while the V band diameter 
showed httle difference from the K band diameter. These results can be interpreted as due 
to the differences of the atmospheric extensions due to the variations of the opacities with 
wavelength. 

More recently, multi-wavelength spatial interferometry with 4 narrow band filters within 
the K band has been done for ji Cep by Perrin et al. (2005). Such an observation is what we 
have been looking for during a long time, and should be regarded as a milestone towards the 
ultimate observations with high resolution both in the spectral and spatial domains. Thus, 
such an observation can be expected to provide a final confirmation of the molecular layers 
outside the photosphere by showing direct evidence for the dependence of the geometrical 
extensions on the molecular opacities. We had to know, however, that the situation is by no 
means so optimistic: The molsphere viewed with the visibility data by Perrin et al. (2005) 
is characterized by the radius of about 1.3 i?* and excitation temperature near 2700 K while 
that viewed with the infrared spectra is characterized by the inner radius of about 2i?* 
and temperature about 1500 K (Tsuji 2000b). Clearly, the discrepancies are too large to be 
attributed to the uncertainties of the observations, and we investigate in this paper if a more 
consistent solution can be obtained from the spectral and visibility data. 

Even for the case of a Orionis, which has been a target of extensive observations with 
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a wide spectral coverage and with a variety of methods, our picture of its outer atmosphere 
appears to have not fully converged yet. An extreme case is the water absorption lines 
observed with a high resolution around 12 /xm by Ryde et al. (2006), who attributed the 
origin of the observed H2O lines to the anomalies of the photospheric structure rather than of 
the outer atmosphere, as already noted before. Such an interpretation is apparently against 
the recent interferometric observations (Weiner et al. 2003; Perrin et al. 2004a) and the 
detailed analyses of the visibility data together with the spectroscopic data (Ohnaka 2004; 
Verhoelst et al. 2006). Such a controversy may be a manifestation of the extreme complexity 
of the outer atmosphere of Betelgeuse, and hopefully be a clue to further progress. 

2. Method of Analyses 

Presently, we have no method of treating the outer atmospheres of cool luminous stars 
consistently. Under such a situation, we introduce an ad hoc model just to provide a frame 
by which numerical analysis such as of the spectra and visibilities can be done. 

2.1. Bcisic Stellar Parameters 

As a boundary condition to the outer atmosphere, we use the classical spherically ex- 
tended LTE model photospheres in radiative and hydrostatic equilibrium. The basic param- 
eters we have applied are summarized in Table 1. We keep most of the basic parameters 
of the photosphere we have used before (Tsuji 2000a,b). It is true that even the effective 
temperatures of red supergiants are not yet well established, and there is no definitive answer 
at present. We assumed Tcs = 3600, K for a Ori, which is consistent with the recent inter- 
ferometry determinations (Dyck et al. 1992, 1996, 1998; Perrin et al. 2004a). We changed 
Tcflf of fi Cep to be 3800, K by the reason to be discussed in Sect. 3. 2. The effect of T^s on the 
predicted infrared spectra can be seen by comparing Fig. 2b and Fig. 3b, for example. 

The CNO abundances are most important, since they give direct effect on the infrared 
spectra. We reanalyzed the equivalent width data of OH and NH measured on the FTS 
spectra of Betelgeuse by Lambert et al. (1984) with our model photosphere of T^s = 3600 K, 
and we confirmed the N and O abundances for the case of T^s = 3600 K determined by 
Lambert et al. (1984), who analyzed the effect of T^s on the derived CNO abundances. As 
for the C abundance, we analyzed some weak lines of the CO first overtone bands measured 
from the high resolution FTS spectra and our C abundance, as shown in Table 1, is a factor 
of two smaller than the value of log Aq/Ah —3.7 (for T^s — 3600 K) by Lambert et al. 
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(1984). The major motivation to have re-analyzed the carbon abundance is the poor fits 
of the CO features at 1.7 and 2.3 fim regions in our previous Figs. 2 & 3 (Tsuji 2000a). 
With our new C abundance, fits in these regions are considerably improved as can be seen in 
Figs.2b-4b, 7b, and 9a-12a. The micro and macro turbulent velocities are also found from 
the same analysis and will be discussed elsewhere with details on the abundance analysis (in 
preparation). The abundance analysis is more difficult for /i Cep, which shows very broad 
lines, and we assume the same abundances as for a Ori. 



2.2. Model Photospheres 

Our model photosphere code is essentially the same as our previous one (Tsuji 1976), 
except that the photosphere is now assumed to be spherically symmetric rather than plane- 
parallel. Also opacity data are updated (see the Appendix of Tsuji (2002a)). We are including 
the radiation pressure Prad a-nd turbulent pressure Ptur in the hydrostatic equilibrium, and 
thus 

Id^ = -g^^, (1) 

p dr 

where p is the density, Pgas is the gas pressure, and 

5'eff = ffgrav " 5'rad " 5'tur (2) 

with 

5'grav — ^ ^ 



G^. (3) 







and 



47r /"^ 

9'rad = — / KuHF^du, (4) 
gtnv = -7—- (5) 

p dr 

The notations have their usual meanings, and we assume Ptur = pCtur^ with the turbulent 
velocity ^tur- The photosphere is stable so long as Qes > and this stabihty limit was defined 
as the "Eddington hmit" for the turbulent plus radiation pressure by De Jager (1984), in 
analogy with the Eddington's well known stability limit for the radiation pressure. 

Usually, integration of model photosphere starts at a very small optical depth in the 
continuum (or optical depth in the mean opacity, e.g. Rosseland mean optical depth tr) such 
as To = 10~^. For T^q — 3600 K and for the other parameters in Table 1, the photosphere 
extends to r(To = 10""^) ^ 730 P© with r(To — 1) ^ 650 P© as shown by the solid lines in 
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Fig.l. However, it is possible to construct a highly extended model photosphere in radiative 
and hydrostatic equilibrium, just starting at a still smaller optical depth within the stability 
limit noted above. An example starting at Tq = 10~^^ (again T^s = 3600 K, and other 
parameters in Tabic 1) is shown by the dashed lines in Fig.l. The photosphere expands to 
r(ro = 10^"'^'^) ~ 1700 i?0 for this case, again with r(ro = 1) ~ 650 i?©. Thus, within the 
framework of the classical model photospheres, a model photosphere can be extended, at 
least formally, to as large as a few stellar radii due to the radiation and turbulent pressures. 
However, it is to be noted that such a model should not be regarded as physically reahstic. 
In fact, the basic assumption such as LTE should certainly be not apphed at such density 
as low as 10~^°grcm~^ (see Fig.l). For this reason, such a model will not be used in our 
actual analysis, but only be used as a reference for an ad-hoc model of the molsphere to be 
discussed in Sect. 2.3. 

If we start the integration from a still smaller optical depth or if we assume a somewhat 
larger turbulent velocity, the photosphere can no longer stay within the De Jagcr's generalized 
Eddington limit and expands without limit. Our purpose here is not to investigate such a 
stability limit, but we hope that such an extended photosphere may give some clues for 
medeling the molsphere. Unfortunately, however, the matter included in the extended part 
(e.g. outside of Tq ~ 10~^) is very small and gives little effect on the spectra as well as on 
the visibilities. It is difficult to deposit more matter in the upper photosphere within the 
framework of the hydrostatic equilibrium model, even though a kind of dynamical effect is 
introduced through the turbulent pressure. One idea may be to consider some kind of shock 
wave by which the rarefied gas can be compressed and heated (e.g. Woitke et al. 1999), and 
the extended photosphere will provide the pre-shock condition for such a model. 



2.3. Ad hoc Model of the Molspheres 

We recognize it difficult to construct a self-consistent model of the molsphere at present, 
and we use an ad hoc model starting from the extended photosphere discussed in Sect. 2. 2. 
We consider a model of molsphere specified with the excitation temperature Tex, gas column 
density A^coi, and inner radius i?in. Then, in the layers above r = i?in in the extended LTE 

photosphere, the temperatures are replaced with Tex and the matter densities are increase 
by a factor iVcoi/A/^LXE) where A'^lte is the LTE column density of the extended photosphere 
above r = The resulting molsphere is isothermal at T^x, with the column density A'^coi, 
and extending between r = i?in and the outer radius r = -Rout, which is near that of the 
starting extended photosphere (e.g. i?out ~ 1700 i?© if starting from the model of Fig.l). 
Since the matter density is rather high near the inner radius, the effect of the outer radius is 
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not important, but outer radius can be changed, if necessary, by starting from the extended 
photospheric model with different Tq ^. The matter densities and temperatures below i?in 
conserve the original LTE values of the extended photosphere, but they give little effect on 
the resulting spectra and visibilities because the matter densities there are very low. 

Our model of the molsphere is not a physical model at all, and no longer in radiative 
and hydrostatic equilibrium, but will be used as a means by which to compute spectra and 
visibilities, and to infer the physical parameters such as temperature, column density and 
size of the molsphere from the observed data. This model is essentially the same as just 
to assume an envelope or shell of given parameters, as is usually done. An advantage of 
our formulation is that the spectral synthesis and related computer codes being used for the 
photospheres can directly be used for our molsphere plus photosphere models with almost 
no changes. 



We compute the specific intensity with our spectral synthesis code for 98 values 

of where /x = cos 9 (9 is the angle between the normal and the direction to the observer). 
Then, the flux iv(r) is given by: 



with the abscissa Hi and weight factor Wi for Gaussian integration ^. Then, F*{r) corrected 
for the dilution effect is obtained by 



where i?* is the stellar radius at r^j ~ 1 and defines the effective temperature T^e through 
L* = AirRlaT^Q (L* is the stellar luminosity and a is the Stefan-Boltzmann constant). 

The spectra are computed with the resolving powers of 5 x 10^ - 10^ and convolved 
with the slit functions matching to the resolutions of the observed spectra. We use the line 



^Notc that i?out differs in the models to be discussed below (Models A - H in Tables 2 & 4) even for 
the same starting extended photosphere, and this is because the hydrostatic equilibrium is solved in each 
starting model with the isothermal molsphere of different Tex- This computations, however, are to determine 
A'lte needed to estimate A^coi rather than to evaluate i?out- 

^We are not using the higher order formula, which does not necessarily result in a higher accuracy, but 
we use the simplest case of n = 2 (e.g. Abramowitz & Stegun 1964) repeatedly. 



2.4. 



Spectra and Visibilities 
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list including H2O (Partridge & Schwenke 1997), OH (Jacquinet-Husson et al. 1999), CO 
(Guelachivili et al. 1983; Chackerian & Tipping 1983), SiO (Lavas et al. 1981; Tipping & 
Chackcrian 1981; Langhoff & Bauschlicher 1993), and CN(Cerny et al. 1978; Bauschlicher 
et al. 1988). All these molecular species (H2O, OH, CO, SiO, CN) are included in the 
calculation of photospheric spectra, but only H2O and CO are included in the molsphere. 

For a spherically symmetric object with the radius /?, the visibility is obtained 
through the Fourier transform of the strip intensity distribution, ^i,{x), defined by 



^,{x)^2 lAx,y)dy, (8) 

Jo 

where Ii,{x,y) is the intensity at {x,y) on the object surface, and the baseline of the inter- 
ferometer is in the x-direction (Michelson & Pease 1921). In the numerical analysis of the 
extended envelope around central star, we follow the formulation we have used before (Tsuji 
1978b): We introduce 

where and d are the angular diameter of the central star and distance to the object, 
respectively, and 

" = -A- 

where I and A are the separation of the two telescopes on the interferometer baseline and 
wavelength, respectively. Then, the monochromatic visibility V^, is given by 

K('y) = / ^u{p)cos{2nvp)dp/ / ^^{p)dp. (11) 
Jo Jo 

For comparison with observed visibilities, which are usually obtained through wide or 
narrow band filters, the monochromatic visibilities are squared and averaged with the filter 
transmissions as weights to have the band averaged visibility, following Perrin et al. (2004a). 



3. The M Supergiant Star Cephei 

With the use of the recent visibility data, we first examine our previous model of the 
molsphere of n Cep based on the infrared spectra alone (Tsuji 2000a,b) and confirm that 
the model is already reasonably consistent with the recent visibility data (Sect. 3.1). We 
further show that the model can be revised to be more consistent with both the visibility 
and spectral data (Sect. 3.2). The resulting model, however, does not agree with the model 
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based primarily on the visibility data (Perrin et al. 2005), and we show that the spectral data 
are as important as the visibility data in modeling the outer atmospheres of cool luminous 
stars (Sect. 3. 3). 



First, we examine a model based on the parameters derived from the infrared spectra 
(Tsuji 2000a,b), as summarized in Table 2 under Model A. It is to be remembered that the 
temperature and column density were estimated from the width and depth of the H2O 1.9 fim 
bands, respectively, on the Stratoscope spectra. The inner radius was estimated from the 
strength of the H2O 6.3 fim bands in emission on the ISO spectrum. Thus, our model consists 
of the photospheric model of Tcs = 3600 K and the molsphere model characterized by the 
excitation temperature Tex = 1500 K, H2O and CO column densities of iVcoi ~ 3.0 x 10^° cm~^, 
and inner radius of i?in ~ 1300 i?© (~ 2R^). The molsphere, starting from the extended 
photosphere of Fig.l, actually extends from R^n ~ 1300 7?© to it!out ~ 1700 i?©, but effective 
contribution comes from the layers close to the inner radius where the density is relatively 
high. 

We now examine this model of the photosphere- moZsp/iere with the recent visibility 
data for the four narrow band regions within the K band by Perrin et al. (2005). We 
compute the intensity lui/j) at the spectral resolution of 0.1 cm~^ for our combined model 
of the photosphere and molsphere, and the strip intensity distributions and visibilities are 
evaluated, as outlined in Sect. 2.4. We apply the angular diameter of the central stellar disk 
of 6^ = 14.11 mas as estimated by Perrin et al. (2005). In evaluating the band averaged 
visibility from the monochromatic visibilities, the filter transmission of each narrow band 
filter is approximated by the Gaussian with the parameters given in Table 1 of Perrin et al. 
(2004b). In comparing the predicted and observed visibihties, value is evaluated by 



where Vi(obs) and (model) are the observed and prdicted visibilities, respectively, Uj is the 
error estimate to V^(obs)^, and N is the data points for each band. 

The predicted visibilities based on the Model A are compared with the observed ones by 
Perrin et al. (2005) in Fig. 2a. The fits are pretty good for the i^239 band which is strongly 
blanketed with the CO and H2O fines {x'^{K239) = 0.23), and for the K215 band relatively 
free from molecular fines {x'^{K215) — 5.67). However, the fits are poor for the K203 band 
which suffers the effect of the H2O 1.9 pim bands (x^(i^203) = 53.68), and for the K222 band 



3.1. Molsphere as Seen by the Infrared Spectra 
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which includes no known strong molecular bands {x'^{K222) — 52.70). The poor fits in these 
bands are largely due to the higher predicted visibilities at the highest spatial frequencies 
compared with the observed ones. Thus, further fine tunings in our modeling should be 
required. 

The predicted near infrared spectrum based on the Model A is compared with the 
Stratoscope spectrum in Fig. 2b. In the following comparisons of the observed and predicted 
spectra, the observed data are shown by the dots and the predicted spectra by the solid 
lines in general. The ordinate scale always refers to the predicted spectrum (in unit of the 
emergent flux from the unit surface area, i.e., in crgcm^^ scc^^ Hz^^), to which the observed 
one is fitted. Also, the predicted spectrum from the photosphere is shown together with that 
from the molsphere plus photosphere in most cases. The details of the Stratoscope spectra 
and the ISO spectra were discussed before (Tsuji 2000a,b), and all the spectra of Cep are 
dereddened with Ay — 1.5 mag. 

This comparison of the Stratoscope spectra with the predicted ones has already been 
done before (Tsuji 2000a), but some improvements are done since then. First, the extra 
molecular layer was approximated by a plane parallel slab in the previous analysis, and this 
was thought to be a reasonable approximation for the near infrared. However, the result 
based on the spherical molsphere shown in Fig. 2b reveals that the sphericity effect is already 
important in the near infrared. This is clearly seen in the CO first overtone bands around 
2.5 fim, whose excess absorption was too large in the previous analysis (sec Fig.3 of Tsuji 
(2000a)) compared with the present result in Fig. 2b. This is because the emission in the 
extended molsphere already has significant contribution to the integrated spectrum. The 
predicted CO absorption around 2.5 /xm, however, still appears to be deeper compared with 
the observed. So far, we assumed that the column density estimated from the H2O features 
applies to CO as well, but it may be possible that the CO column density is somewhat smaller. 
But we assume that the column densities for CO and H2O are the same for simplicity. 

Second, the region of the CO second overtone bands around 1.7 /xm, which was depressed 
too much compared with the observation in our previous analysis (Fig.3 of Tsuji (2000a)), 
is now improved in Fig.2b, and this is the effect of the revised carbon abundance noted 
in Sect. 2.1. Also, the improved fit near 2.5 /im is partly due to the revision of the carbon 
abundance as well. As a result, the near infrared spectrum observed by Stratoscope including 
the H2O 1.4 and 1.9 fim bands can be well matched with our prediction based on the Model 
A. Although the observed fiux short-ward of 1.3 /im shows some excess compared to the 
prediction, this may be a problem of the photospheric spectrum, which is also shown in 
Fig. 2b and shows only CO and CN bands. 

The predicted spectrum in the H2O 1/2 band region based on the Model A is compared 
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with the ISO spectrum in Fig.2c. This comparison has also been done before with the use 
of the spherically extended molsphere model ^ (Tsuji 2000b), but we now use the H2O line 
list by Partridge & Schwenke (1997) instead of the HITEMP database (Rothman 1997) used 
before. The tendency that the H2O 1^2 bands appear in emission can be well reproduced 
with our Model A, but the predicted emission is still not strong enough to explain the 
observed spectrum in the region around 6.8 /xm. The predicted spectrum based on the 
model photosphere appears at the lower left and the most hnes seen are due to the tail of 
the CO fundamental bands. 



3.2. Some Fine Tunings and a Resulting Model 

Inspection of Figs. 2a-c reveals that the predicted visibilities and near infrared spectra 
based on the Model A are not yet fully consistent with the observations. The major problem 
is that the predicted visibilities at the K203, especially at the highest spatial frequency, are 
much higher than the observed. Also, the predicted emission of the H2O1/2 bands is not 
sufficiently strong to be consistent with the observation. Such inconsistencies can be relaxed 
by changing (a) the excitation temperature of the molsphere, (b) the inner radius, and/or 
(c) the gas column density. However, the column density is already well consistent with the 
near infrared spectrum and we will not consider the possibility (c) further. Also, since the 
visibilities based on the model A are already well consistent with the observations at K215 
and i^239 bands, the possibilities (a) and (b) must be considered so that the good fits in 
the visibilities already achieved can be conserved. We examine the possibilities (a) and (b) 
systematically by changing the temperature by a step of lOOK and the inner radius by a step 
of 0.1 R^. Also, as to the observed excess in the short- ward of 1.3 //m, we notice that the 
recent angular diameter measurements suggested a higher effective temperature of 3789K 
(Perrin et al. 2005), and we apply Teg = 3800 K in all the models to be discussed below, 
instead of 3600 K used in the Model A. The extended phtosphere of T^s = 3800 K, used as 
the starting model for the molsphere, is quite similar to the case of T^e = 3600 K shown in 
Fig.l. 

First, we examine the effect of changing i?in from 2.0 i?* ~ 1300 i?0 of the Model A, 
keeping Tex to be 1500 K. After some trials and errors, we find that a case of Rm ~ 2.2 i?* ^ 
1430 i?0 shows some improvements. First, the resulting predicted visibilities for the four 
narrow band regions are compared with the observed data in Fig. 3a: The fits for the A'203 
and /\ 222 bands arc somcw hat improved with x^{K203) = 40.51 and x^{K222) = 40.95, but 



''This model was a preliminary version of the Model A and not exactly the same with it. 
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still not very good. On the other hand, fits for the 7^215 and ^239 bands remain reasonable 
with x'^{K215) = 4.66 and x'^{K239) = 13.61. The good fits in the visibilities for the i^239 
band with the Model A is somewhat degraded but still remain to be reasonable. Second, 
the predicted flux in the region short-ward of 1.3 ;um is increased as a result of changing 
Tcff of the photosphere to 3800K from 3600K, and the fit in this region now appears to be 
reasonable as shown in Fig.3b. Finally, the observed and predicted spectra in the 6 - 7/xm 
region are compared in Fig. 3c. The predicted emission around 6.8 /im increases appreciably 
and it is now sufficiently large to account for the observed emission as shown in Fig.3c. We 
refer to this model as Model B and its major parameters are summarized in Table 2. 

Although the model B is considerably improved compared to the model A, we examine 
another possibility of changing Tgx of the molsphere from 1500 K of Model A while i?in remains 
to be 2.0 i?*. The resulting predicted visibilities for a case of Tex = 1600 K are again compared 
with the observed ones in Fig.4a: The fits for the K203 and K222 bands are both improved 
compared to the results for the model B, with x^(^203) = 29.24 and x^(^222) = 27.76. At 
the same time, fits for the K239 band is further degraded but still remain to be acceptable 
with x^(^239) = 17.68. The fits for the K215 band remain fine with x^{K215) = 4.52, 
and this result that the fits for the 7^215 band remain nearly unchanged may be a natural 
consequence that this band region is a good continuum window. We conclude that this model, 
to be referred to as Model C, provides better fits to the observed visibilities compared to 
the Model B. The predicted near infrared spectrum shown in Fig.4b shows almost no change 
from Fig.3b, and we confirm that the 1.4 and 1.9 //m absorption bands depend primarily on 
the column density which is unchanged throughout. The predicted emission around 6.8 /im 
reasonably account for the observed emission as shown in Fig.4c. The major parameters of 
the Model C are summarized in Table 2. 

We examined several other models around the models B and C, but no significant 
improvements were obtained. For example, a case of increasing i?in to 2.1 R^, keeping Tgx — 
1600 K, results in degrading the fits in visibilities, although fits in spectra remain nearly the 
same. A case of increasing Tex to 1700 K, keeping i?in = 2.0 i?*, also results in degrading the 
fits in visibilities, while fits in spectra remain nearly the same. If -Rin is decreased to 1.9 7?* 
at Tgx = 1700 K, the fits in visibilities for the K222 band are somewhat improved*^ compared 
to those for the Model C, but the predicted emission around 6.8 /xm is reduced to the level 
of the Model A. On the other hand, if Tex is decrease to 1400 K, R^^ had to be increased to 
2.2 R^ to keep the reasonable fits for the spectra, but the fits for the visivihties are degraded 
appreciably. 



^The values for these fits are: x^{K2Q2,) = 35.64, x^{K2lb) = 8.91, x^{K222) = 18.04, and 
X^(ii'239) = 31.15 
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We conclude that the Model C is a possible best model that reasonably accounts for 
both the observed visibilities and infrared spectra simultaneously, within the framework 
of our highly simplified molsphere model. It is encouraging that our model could repro- 
duce the general tendency that the visibilities in the region with strong H2O and/or CO 
bands {K203, K239) are lower than those in the region with less molecular absorption 
{K215, K222), in agreement with the observations, even though the fits are not very good 
for the K203 and K222 bands. This result imphes that the molsphere model based on the 
infrared spectra alone is reasonably consistent with the visibility data not known at the time 
when the modeling was done. This may be because the infrared spectra already include 
some information on the spatial structure of the extended atmosphere. It is most important, 
however, that the model is now examined directly with the interferometry, which is a direct 
probe of the spatial structure of the astronomical objects. 



3.3. Molsphere as seen by the Visibility Data 

We notice, however, that the model parameters of our Model C show significant differ- 
ences with those derived from the visibility data by Perrin et al. (2005) themselves: Their 
excitation temperature is near 2700 K compared with 1600 K of our Model C and their radius 
of the molecular layer is about 1.3 i?^ compared with the inner radius 2.0 -R^, of our Model C. 
To investigate the origin of the differences, we think it useful to analyze their model in the 
same way as in our models. For this purpose, we refer to a model based on their parameters 
as Model D. Since their molecular layer is assumed to be a thin shell, we design the Model D 
in our formulation to be extending from i?in = 851 to it!out — 858 i?*, for which the initial 
extended photosphere starts at Tq = 10~^, and thus the effective location of the thin shell of 
the Model D is at about r 1.3 R^. 

The gas column densities were not provided in the model of Perrin et al. (2005) but a 
mean optical thickness for each narrow band region was given. Accordingly, we transform 
their optical thickness at each filter band region to the column densities with the absorption 
cross sections of CO and H2O evaluated from the line hst mentioned in Sect. 2. 4 and shown 
in Fig. 5, where the high resolution results are given with the black while the smeared out 
straight means by the white lines. The resulting column densities summarized in Table 3 
appear to be different for different bands, while this should be the same at least for 2.03, 
2.15 and 2.22 /im in which H2O is the major source of opacity. The column density of CO 
can of course be different from that of H2O. We assume that the column densities of CO 
and H2O for the 2.39 //m region are the same as in Table 3, just as an example. But the 
optical thickness of 3.93 for the 2.39 //m region given by Perrin et al. (2005) is large enough to 
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completely mask the photospheric spectra and the column density corresponding to r = 3.93 
should anyhow be very large. 

The major parameters of the Model D are summarized in Table 2 and we can now 
proceed as in our Models A - C with these parameters. The resulting visibilities are shown 
in Fig. 6 together with the observed data by Perrin et al. (2005). The predicted visibilities 

using the column densities as input parameters nearly reproduce the ones for the K215, K222, 
and K239 bands shown in Fig.l of Perrin et al. (2005) using the effective optical thicknesses 
as input parameters. However, such an agreement cannot be found for the K203 band, and 
this means that the transformation of the mean optical depth to the column density could 
not be done correctly. We cannot understand the reason for this, but the H2O cross-section 
at temperature as high as 2700 K may not necessarily be accurate enough for the region 
of the K203 band. Thus, we may conclude that the visibility computed with the mean 
optical thickness can be a reasonable approximation to the one based on the monochromatic 
visibilities evaluated with the realistic line list at high spectral resolution, only if the optical 
thickness correctly reflects the detailed line opacities. ^. 

As shown in Fig. 6, the predicted visibilities based on the Model D show fine fits with the 
observed ones for the K215, K222, and K239 bands, for which x^(K215) = 4.02, x^(X222) = 
14.23, and x^(-^239) = 2.32. This is what can be expected, since our visibilities agree well 
with those by Perrin et al. (2005), which are already known to be well consistent with 
the observed data. For the same reason, the fits for the A'203 band are quite poor with 
x'^{K203) = 47.93. It is, however, possible to improve the fits in the K203 band by just 
changing the column density somewhat, and we find a reasonable fit with x^(i^203) = 16.99 
for A^coi(H20)= 4.2 x 10^°cm~^ after a few trials and errors (see dotted line in Fig. 6). If 
we can choose four values for the free parameter such as the mean optical depth or the 
column density for the four observed data, it is certainly possible to have good fits for all 
the observed data. However, such fits cannot be regarded as justification of such a model as 
D, since the values of the column densities are so different for the different bands (see Table 
3). 

Since the column densities should be unique for the different bands, we modify the 
Model D so that A^coi(H20)= 2.8 x 10^° cm'^, just as an example, for all the 4 narrow band 
regions. Then, to account for r = 3.93 of the K239 band which includes CO and H2O, it 
turns out that A^coi(CO)= 9.3 x 10^^ cm~^ with the absorption cross-sections given in Table 



^Furthermore, this result can be appUcd to the case that the straight mean opacity is a reasonable 
approximation as for the hot water vapor. Note that this result can no longer be applied to CO for which 
the straight mean opacity no longer describes the spectrum accurately 
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3. For this modified Model D, to be referred to as Model D*, we first evaluate the visibilities 

and the results arc again compared with the observed data in Fig. 7a. The fits are generally 
fair with x^(X203) = 17.78, x^(^215) = 25.24, x^(^222) = 14.21, and x'^{K239) = 30.45. 
The fits for the K203 and K222 bands are better than those for the Model C while the fits 
for the K215 and 7^^239 bands are better in the Model C. Thus, it is difficult to decide which 
of the Model C or Model D* is to be preferred from the visibility analysis alone. ^ 

An advantage of using the column density instead of the mean optical thickness is that 
the spectra can be evaluated for the same input parameters. The near infrared spectrum for 
the Model D* is compared with the Stratoscope data in Fig. 7b. The fit is rather poor and 
the other A^coi(H20) values in Table 3 produce too weak or too strong H2O absorption bands. 
We also evaluate the spectrum in the 6 - 7 iim region for iVcoi(H20)= 2.8 x 10^°cm~^ and 
compared with the ISO data in Fig. 7c. The predicted spectrum still appears in absorption 
and cannot be matched with the observed one at all. The results are more or less the same 
for other values of A^coi(H20) in Table 3. 

In our fine tunings in Sect. 3. 2, we increased the temperature from 1600 to 1700 K and, 
at the same time, decreased the inner radius to from 2.0 to 1.9 i?*, and obtained reasonable 
fits to the observed visibilities. If we pursued a solution in this direction, we might increase 
the temperature further and decrease the inner radius at the same time. The resulting 
model might be similar to the Model D*, but we might reject such a case since such a high 
temperature - small size model might not explain the 6 - 7/im emission as can be inferred 
from the result shown in Fig. 7c. Perrin et al. (2005) did not consider the spectral data 
and might be led to the high temperature - small size model that might have satisfied the 
visibility data well. 

We conclude that it is difficult to decide the unique solution for the molsphere model 



^Compared to the Model C, the fits between the predicted and observed visibihties are improved for the 
i('203 and K222 bands in the Model D*. This is because the tails of the H2O 1.9 and 2.7/um bands are well 
excited in the regions of the and K222 bands, respectively, at a rather high Tgx of 2700 K (Fig. 5), and 

the visivilities in these bands arc reduced to be better agreement with the observed visibilities. However, 
the same effect makes the region of the K21^ band too opaque and the predicted visibilities are reduced too 
much to be fitted well with the observed data. For this reason, the fits for the K21f) band are wrose in the 
Model D* tnan in the Models B and C. For the K2?>'\} band, the fits arc also wrose in the Model D* tnan in 
the Models B and C. Here, A^coi(CO) is chosen to be consistent with r — 3.93 for the K2'i% band region. A 
problem, however, is that the spectrum of CO, unlike the case of H2O, is not well smeared out and hence 
does not block the radiation so effectively with the column density of A^coi(CO)= 9.3 x 10^^ cm~^ derived via 
the straight mean opacity. Remember that the case of A^coi(H20) = A^coi(CO)= 3.3 x 10^^ cm~^ assumed 
in Fig. 6 could predict a lower visibility, and this is because the H2O spectrum is well smeared out and the 
column density was estimated properly via the straight mean opacity. 
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from the visibility data alone, and it is indispensable to apply the spectroscopic data at the 
same time to have some idea on the model of the stellar outer atmosphere. At the same time, 
it is difficult to have good fits to the observed visibilities at all the spatial frequencies within 
the framework of the simplified uniform spherically symmetric models of the molsphere. 
Certainly, some details of the fine structure of the outer envelope should be considered to 
have better fits at higher spatial frequencies, and our present analysis is only a very initial 
stage on the interpretation of the visibility data. 

4. The Case of a Orionis 

Bctelgeuse is an object probably best observed with a variety of methods and new 
observed data are still accumulating. At the same time, Bctelgeuse is a quite complicated 
object, and how to interpret the observed data is still controversial at least partly. We will 
show that unexpected water lines observed in Bctelgeuse may be difficult to be explained 
by anomalies within the framework of the photospheric models (Sect. 4.1). On the other 
hand, the infrared spectra and visibility data consistently show the presence of the extra 
molecular layers and provide consistent estimations of their basic parameters (Sect. 4.2). In 
the mid infrared region, some recent observations subjected to controversy, but we look for 
a possibility to relax such an issue (Sect .4.3). 

4.1. Can a Cool Photosphere Explain the Infrared Spectra of Betelgeuse? 

This is a question addressed recently by Ryde et al. (2006), who observed high resolu- 
tion spectra of Betelgeuse in the 12 //m region and suggested that the observed H2O lines 
can be interpreted with the anomalous structure of the photosphere rather than assuming 
the presence of the extra molecular layers. They showed that their high resolution 12 fim 
spectra could be well fitted with the predicted ones based on the cool photospheric models. 
However, it can be shown below that the near infrared spectra cannot easily be matched by 
an anomalous structure of the photosphere. 

Ryde et al. (2006) suggested that the observed H2O pure rotational lines in Betelgeuse 
can be accounted for with the thermal structure approximated by the classical model photo- 
sphere of Teff ~ 3250 K and that their observations can be explained with a rather small H2O 
column density of A^coi ~ 5 x lO^^cm"^. This is possible because the /-values of the H2O 
pure rotational fines are pretty large. In Fig. 8a, the spectra of the H2O pure rotation fines 
predicted from the classical spherically extended model photospheres of — 2800, 3000, 
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3200, 3400, and 3600 K (M = ISM©, R = 650 i?©, ^^icro = 5kms-^ ^^acro = 10kms-\ 
throughout) are shown. Thus, weak H2O pure rotational hues can be observed already at 
Teflf = 3600 K and increasingly stronger towards the lower T^a 's. At Teg ~ 3200 K, H2O lines 
can certainly be observed even if they are diluted by the silicate dust emission in Betelgeuse. 

However, with such a small column density, it is not possible to explain the H2O 1.4 
and 1.9 /xm bands observed with Stratoscope II, since the /-values of the near infrared H2O 
lines are 1-2 orders of magnitude smaller than those of the pure-rotation lines observed 
in the 12 /xm region. In Fig.Sb, we show the near infrared spectra predicted from the same 
models as used for Fig. 8a and reduced to the resolution of the Stratoscope observations. It 
can be confirmed that the II2O 1.4 and 1.9 /xm bands cannot be seen in the models with 
Teflf = 3200 K and higher, but appear in the models of T^q — 3000 K and lower. It is clear 
that the H2O 1.4 and 1.9 //m bands seen on the Stratoscope spectra (Figs. 9a - 12a) cannot be 
accounted for with the cool photosphere of Teflf ~ 3250 K. Ryde et al. (2006) argued that the 
photosphere of Betelgeuse can be very peculiar, but the relative behaviors of the near and 
mid infrared spectra should remain nearly the same so long as any peculiarity remains within 
the photosphere. We conclude that the near and mid infrared spectra cannot be consistently 
interpreted by a cool atmosphere wihin the framework of the classical model photosphere. On 
the other hand, Ryde et al. (2006) suggested a possibility of an inhomogeneous atmosphere 
with a cool component giving rise to the 12 /xm water lines and a hot component giving rise 
to the near-infrared H2O bands. Such a possibility should hopefully be examined further 
with other arguments suggesting the inhomogeneity in red supergiant stars (Sect. 5.1). 



4.2. Can a Molsphere Explain the Infrcired Spectra of Betelgeuse? 

For Betelgeuse, we have assumed the presence of an extra molecular layer above the 
photosphere, and estimated Tgx = 1500 ± 500 K and A^coi = 10^° cm~^ from the Stratoscope 
spectra based on a simple slab model (Tsuji 2000a). Also, we have suggested that the weak 
6.6 //m absorption observed with ISO can be explained if the inner radius of the molsphere 
is it!in ~ 1.5ii!* (Tsuji 2000b). This result, however, was given without a detailed numerical 
analysis and may not be accurate enough. We now re-examine such a model of Betelgeuse in 
some detail with some improvements such as noted in the analysis of fi Cep. The observed 
spectra of Betelgeuse are discussed in detail before (Tsuji 2000a,b) and dereddened with 
Ay = 0.5 mag. 

For this purpose, we recall that the near infrared spectra observed with Stratoscope II 
are relatively free from the effect of spherical extension of the molsphere, and thus they are 
very useful probes of the column density and temperature. Especially, we believe that the 



-18- 



column density of H2O can be well fixed from the strengths of the 1.4 and 1.9//m bands to 

be iVcoi ~ 10^°cm~^. We confirm that this result remains unchanged for possible values of 
Tex and -Rin in the spherically extended molsphere models. Then, although the H2O 6.3 /xm 
bands do not appear in emission but only in weak absorption in the case of a Ori, they can 
be used to constrain the extension of molsphere and to refine the temperature. Thus the 
nature of the extra molecular layer can also be inferred from these infrared spectra alone as 
in the case of fj, Ccp. The result should further be examined in the light of the interferometric 
observations. For this purpose, we refer to the recent work by Perrin et al. (2004a), who 
suggested a 2055 K layer located at 0.33 i?* above the photosphere from the visibihty data 
at the K, L, and 11.15 bands. 

We assume T^s — 3600 K for the central star and start from the extended photosphere 
of Fig.l. We examine first the case of Tex = 1500 K and systematically changed the inner 
radius i?in of the molsphere from 1.3 R*, approximately the radius of the molecular layer by 
Perrin et al. (2004a). The near infrared spectrum can be fitted reasonably well for any values 
of Rin up to ^ 2/?* with Nco\ = lO^'^cm"^, but the 6 - 7/xm spectrum appears in strong 
absorption for i?in ~ 1.3 i?* and weakens towards larger i?in. It finally turns to emission 
for ~ 1.9 -R*. After all, reasonable fits in the 6 - 7/im region as well as in the near 
infrared spectra can be obtained for i?in ~ 1.7 i?* (Model E) as shown in Fig. 9. However, 
the predicted spectrum shows appreciable excess in the region long- ward of 7 /im compared 
with the ISO spectrum as shown in Fig.9b. 

Next, we examine the case of T^x = 1750 K, and reasonable fits in the near infrared 
and 6 - 7/i.m regions are obtained for i?in ~ 1.5 i?* (Model F) as shown in Figs. 10a and b, 
respectively. This case is close to the Verhoelst et al. (2006) 's model except for the column 
density. We examined their value of A^coi = 2x10^^ cm~^ and the resulting near IR spectrum 
is shown by the dotted fine in Fig. 10a. This column density may be too small to explain the 
H2O 1.4 and 1.9/xm bands. The dip at 6.6 //m due to the H2O 1^2 bands can be reasonably 
accounted for but the disagreement in the region long- ward of 7 /im cannot be resolved as 
shown in Fig. 10b. 

For the case of Tex = 2000 K, a possible best fits are obtained for R^n ~ 1.4 i?* (Model 
G) as shown in Fig.ll. This case is close to the Ohnaka (2004)'s model except for the column 
density, and we examined his value of iVcoi = 2 x 10^° cm~^. The resulting near IR spectrum 
is shown by the dotted line in Fig. 11a, and this column density may be too large to account 
for the H2O 1.4 and 1.9 /xm bands. The near infrared spectrum is relatively insensitive to 
Tex and i?in, but depends critically on iVcoi. The predicted excess in the region long- ward of 
7 iim compared with the ISO spectrum is somewhat relaxed but the excess still remains as 
shown in Fig. lib. 



-19- 



With a hope to improve the fits further, we examined the case of Tex = 2250 K. In this 

case, ^ 1.3 -R* (Model H) results in reasonable fits as shown in Fig. 12. Although the fits 
do not show any drastic improvements, inspection of the cases of Tgx from 1500 K to 2250 K 
reveals that the fits in some details (e.g. strengths of the absorption features, overall shape 
of the spectra etc.) are generally better for the higher Tex, as can be seen through Fig. 9 to 
Fig. 12. If we further increase Tex to 2500 K, however, absorption and emission in the 6-7 /^m 
region just cancel for ii!in ~ 1.3 i?*. Thus, we may stop our survey here, and we summarize 
the major parameters of the 4 models discussed above in Table 4. 

For a Orionis, the visibility data at the narrow bands as for ^ Cephei are not yet known, 
but those at the wide band K filter are available (Perrin et al. 2004a). We evaluate the band 
averaged visibilities for the standard K band filter as for the narrow band filters discussed 
in Sect. 3, and we assume that the angular diameter of the central stellar disk is — 42 mas 
following Perrin et al. (2004a). First, we compare the visibihty curve predicted for the model 
photosphere of Tcfj = 3600 K (tq = 10~^) with the observed visibility data at the K band in 
Fig. 13, and find that it does not fit well to the observed data (x^ = 300.99). The case for 
the extended photosphere (tq = 10^^"^) discussed in Sect. 2. 2 shows only minor changes, and 
the extended as well as the classical photospheres cannot explain the observed visibilities at 
the K band. 

Next, we compare in Fig. 13 the predicted visibilities at the fT-band for the 4 models 
(Models E - H in Table 4) discussed above, and the values are 134.26, 102.21, 29.56, and 
16.18 for the Models E, F, G, and H, respectively. The fits in the first lobe of the visibility 
curves are better for the models with the higher Tex and smaller R\a (i.e. Model G and 
especially H in Table 4). However, the fits in the second and third lobes appear to be more 
consistent with the models of the lower Tex and larger R^^ (i.e. Models E and F in Table 
4). The first lobe may reflect the basic structure of the object while the second and third 
lobes may depend on some fine structures. Thus, considering the consistent results from 
the infrared spectra and visibilities, we may conclude that the molsphere of Betelgeuse is 
characterized by the rather compact size {R-m ~ 1.3 i?* ), rather high temperature (Tex ~ 
2250 K) and modest column density (A^"coi ~ 10^°cm~^) (i.e. the Model H of Table 4). 
Anyhow, it is clear that the visibility data can be interpreted more consistently with the 
molsphere around the photosphere. We conclude that a molsphere explains not only the 
infrared spectra (the mid IR spectra are discussed in Sect. 4.3) but also the K band visibihty 
data consistently. 
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4.3. What the Mid Infrared Data Tell Us? 

The region short- ward of 7.5 /im discussed in Sect. 4. 2 may be almost free from the effect 
of dust, which, however, will have significant effect in the longer wavelength region. First, we 
examine the SiO feature near 8 /xm, which was used as a further check of the molecular layers 
by Verhoelst et al. (2006). The band-head region of the SiO fundamentals observed by ISO 
can be reasonably fitted by our prediction with A^coi = 10^° cm^^ as shown in Fig. 14, and this 
column density may not be high enough to mask the photospheric spectrum. However, the 
fit in the long- ward of 7.6 //m is quite poor, and the ISO spectrum shows a large depression 
centered at about 8 /im compared with our prediction. This depression may partly be due to 
the effect of SiO in the molsphere, which we have not considering yet. For simplicity, we have 
considered only CO and H2O in the molspherehut the effect of other molecules including SiO 
should certainly be examined. The effect of dust should still be minor in Fig. 14, except for 
a possible tail of the alumina emission suggested by Verhoelst et al. (2006) in the long-ward 
of 8 /im. 

In the 10 /im region, some interesting observations both on spectra and visibilities have 
been done recently. The interferometric observation in the 11 /xm region revealed that the 
resulting apparent diameter of a Ori is about 30% larger than those measured in the near 
infrared (Weiner et al. 2003). This observation has been done in narrow wavelength bands 
which are apparently free from molecular lines, and high resolution spectra confirmed that 
there is no significant spectral feature in the bands. Weiner et al. (2003) interpreted that 
this result should be due to the effect of the continuum opacity, and that a possible presence 
of hot spots may be responsible to the reduction of the near-IR apparent size. On the other 
hand, Ohnaka (2004) showed that the apparent large size can be interpreted as due to the 
H2O opacity of the extended envelope while H2O absorption lines appear to be weakened 
because of the filling in by the emission of the extended envelope. In view of the effect of 
molecular bands in the K band region in a Ori as well as in Cep, this interpretation may 
be reasonable so far as the H2O layer is concerned. Also, the effect of silicate dust shell 
located far from the stellar surface on the spectra as well on the visibilities has fully been 
taken into account. 

Recently, Ryde et al. (2006) observed the 12 //m H2O fines with a higher resolution 
compared to the previous observations by Jennings & Sada (1998) and argued that the 
lines should be originating in the photosphere rather than in the extra molecular layers. As 
a support for this, Ryde et al. (2006) showed that the published realizations of molspheres 
cannot reproduce the observed high resolution spectra in the 12 /xm region but rather predict 
the 12 fira H2O lines in emission. But whether the lines appear in emission or in absorption 
depends on the parameters assumed, and absorption spectrum also appears as shown by 
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Ohnaka (2004). Further test, however, can be done by the details of the spectra such as 
the hnc profiles and relative intensities, which may differ for optically thin case suggested 
by Ryde et al. (2006) and for optically thick case of the spherically extended molecular 
layers: Ryde et al. (2006) showed that the observed details of their high resolution 12 fim 
water spectra fit well with their optically thin predictions, and this result will add a strong 
constraint on the model of the outer atmosphere as will be discuused below. 

In this connection, we may examine the effect of a new dust shell suggested by Verhoelst 
et al. (2006), who noticed that there is an additional component in the ISO spectrum in the 
region extending from 8 to 17 fira with a peak at 13 fim and identified it with the emission of 
the amorphous alumina (AI2O3). Verhoelst et al. (2006) showed that the 11 /xm visibiUtics 
can be better fitted with this new dust component in addition to the H2O layer. The alumina 
shell and the H2O layer are located at about the same height above the photosphere according 
to Verhoelst et al. (2006), and it is quite possible that only a part of H2O in the molsphere 
will form an optically thin layer above the alumina shell, which may act as a continuum 
background source for H2O to form absorption. Thus, the observed 12 ^m absorption lines 
may indeed be produced in an optically thin layer as suggested by Ryde et al. (2006), but 
this does not imply that they are originating in the photosphere. 

Summarizing, water layer itself may have a rather large column density of about 10^° cm"' 
as almost uniquely determined from the Stratoscope spectrum, which is also almost free from 
the effect of dust. But the mid infrared water lines may be produced in an optically thin 
layer which is only a part of the whole water layer. This explains not only the argument 
of Ryde et al. (2006) that the mid-infrared water lines are formed in an optically thin layer 
but also the rather small H2O column density of 2 x lO^^cm"^ determined from the 11 //m 
spectrum by Verhoelst et al. (2006). The column density of 5 x lO^^cm"^ by Ryde et al. 
(2006), however, is appreciably smaller than that by Verhoelst et al. (2006). This can be 
explained by the effect of the emission from the extended part of the molecular layer as 
detailed by Ohnaka (2004): Verhoelst et al. (2006) might have considered this effect in their 
extended model of r ^ 1.4 i?^, and then a larger column density is needed to correct for the 
weakening due to the emission component, while Ryde et al. (2006) assumed photospheric 
origin without such an effect and then a smaller column density is sufficient just to account 
for the observed absorption. 

As noted by Verhoelst et al. (2006) themselves, however, the presence of the amorphous 
alumina cloud is still a hypothesis that consistently account for the spectral and visibility 
data, and further confirmation should certainly be required. In recognizing dust in astro- 
nomical objects, it is a general difficulty that dust shows no clear spectral signature, but 
interferometry could recognize the dust cloud with no spectral signature through its geomet- 
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rical extension. An interesting example is the A^-band spectro-interferometric observation of 
the sihcate carbon star IRAS 08002-3803 by Ohnaka et al. (2006), who suggested the presence 
of a second grain species in addition to sihcate around this sihcate carbon star. Thus further 
detailed interferometry in the mid-infrared region will be quite useful to clarify the nature of 
the dust species around Betelgeuse. Also, the suggested temperature of 1900 K is too high 
for AI2O3 to condense in thermal equilibrium at low densities, and some non-equilibrium 
processes may be required. It is to be noted that even the formation of H2O has the same 
difficulty of how it is possible to form water in the rarefied outer atmosphere. 

In conclusion, by virtue of the alumina dust shell recently suggested on the ISO spectra 
by Vcrhoclst ct al. (2006), we found a posibility that all the mid infrared observations 
including the spectra and visibility can be interpreted consistently at last, but this is possible 
only if we assume the presence of the extra components, consist not only of the molecular 
layer but also of the dust shell beyond the photosphere. Our present model consisting of 
gaseous molecules alone is certainly not applicable to the mid infrared region and, for this 
reason, we confine o\ii analysis to the region below 7.5 //.m where is almost free of the effect 
of dust. The next step should certainly be to include warm dust such as alumina as an 
important ingredient in the outer atmosphere, and we hope that more observations and 
theoretical analyses for this purpose could be developed. 



5. Discussion 
5.1. Modelings 

Given a model, any observable such as spectrum and visibility can be calculated directly, 
but the reverse problem to specify a model uniquely from the observed data is not necessarily 
so straight forward in general. This is because an observed result generally depends on 
several model parameters. We could determine a set of parameters that are consistent with 
the known observed data, but this is by no means a unique solution and may be regarded 
as a possible one solution at best. It is to be remembered that even the definitions of the 
parameters are different by the different authors. 

In fact, the definition of the parameter that specifies the size of the molsphere are 
different. For example, a thin shell model with no geometrical thickness is characterized by 
a single parameter - the radius of the thin shell (e.g. Pcrrin ct al. 2005). If the molecular 
layer extends from the stellar surface to a certain height, a single parameter that specify the 
outer radius is sufficient (e.g. Ohnaka 2004). The shell model of a finite thickness requires the 
inner and outer radii (e.g. Verhoelst et al. 2006). Our models emphasize the inner radius. 
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where the density is the largest, reflecting the starting model photosphere in hydrostatic 
equilibrium. Nevertheless the outer radius i?out may have some effect, but we have not 
examined such an effect in detail in this paper However, we think it not useful to increase 
the number of parameters that may not be very essential. 

Detailed parameter fittings may not be our final purpose, but it may be of some interest 
to know if the molecular layer is very thin or rather thick, or if it is detached or not from the 
photosphere. Such information can be of some help in considering the origin of the molecular 
layer, but it may be difficult to decide such details of the geometrical configurations from 
the interferomctric observations at present. However, combined with the spectroscopic data, 
some information can be obtained. One interesting possibility for this purpose is the case 
that the gaseous molecules and dust cloud coexist as in the case of water and alumina clouds 
suggested by Verhoelst et al. (2006). In such a case, the column densities of water estimated 
from the spectral regions with and without dust background may give a clue for the extent 
of the water layer (cf. Sect. 4.3). 

As already discussed in Sect. 3.3, it was not possible to have good fits to the observed 
visibilities at all the spatial frequencies, and this fact implies that our present modeling is 
too simplified. Also, our modeling is still a kind of parameter fitting for this simplified model 

and not a physical modeling yet. Thus we can have a rough idea on the outer atmosphere 
of red supcrgiant stars, but this is the present limitations in our modelings. Certainly, our 
modelings of the available data is to find some chics for a more physical modeling, but it 
is only recently that the visibility data are made available and we arc just starting towards 
such a purpose. Certainly, more interferometry data for a wider coverage in spectral region 
as well as in spacial frequency are highly important. 

We assumed the presence of a rather warm gaseous cloud above the photosphere, and 
the resulting model of m,olsphere may be rather massive. For example, the H2O column 
density of 3 x 10^° cm^^ estimated for fi Cep may imply the hydrogen column density of the 
order of lO^^cm"^, if the oxygen abundance in Table 1 can be assumed, and the total mass 
of the molsphere of /i Cep is as large as 10"^ Mq. This will be sufficient as a reservoir for the 
gaseous mass-loss outfiow of about lO^^Mgyear"^ (Josselin et al. 2000). Also, the hydrogen 
envelope can be large enough for the radio continuum to be opaque, and such an effect is 
already observed in the radio domain of Betelgeuse (Lim et al. 1998). 



^For example, if i?out is decreased in the Model C for ^ Cep, then Rm must be increased a bit and/or 
Tex must be changed slightly, to maintain the fits in the visibilities as well as in the 6.6 /xm emission. Thus 
many solutions can be possible around a solution such as the Model C. At present, however, there may be 
no means by which to discriminate such minor difference in these different solutions. 
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Now a more serious problem is how such a new molecular layer can be formed and how 
it can be stable in the outer atmosphere. We propose it as an observational requirement and 
physical interpretation had to be deferred to future studies. However, outer atmospheres of 
red supergiant as well as giant stars are not yet well understood and it is no wonder that 
such a new feature had to be introduced. In fact, more or less similar features are already 
known: For example, it has been well accepted that the high temperature gaseous envelope 
referred to as chromosphere exists in the outer atmosphere of cool stars. However, how the 
chromosphere is formed and how the gaseous matters are transferred to the higher layers 
above the photosphere have not been known for a long time. In fact, it is possible that the 
molsphere has a close connection to the chromosphere, and that they may simply represent 
cooler and hotter phases of the same phenomenon. From this point of view, a possible outer 
atmospheric origin of water detected in the K giant star Arcturus by Ryde et al. (2002) may 
not necessarily be excluded, since the chromosphere is known in this K giant star. Also, 
in cooler supergiant stars known as maser sources, the presence of huge water clouds is 
known as an observational fact (Sect. 5.3). Thus, we think it useful to accept the presence 
of molspher and to investigate its nature with all the available observational and theoretical 
possibilities. 

Finally, we must also remember that even the photosphere of red supergiants includes 
many unsolved problems. For example, some infrared lines of OH in Betelgeuse show anoma- 
lous intensities that cannot be explained by the photospheric models (e.g. Lambert et al. 
1984), and the origin of the super-sonic turbulent velocities is also unknown. One interesting 
problem often mentioned is an inhomogeneity in the photosphere. A possibility of a large in- 
homogeneity was suggested by Schwarzschild (1975), who argued that only a few convective 
cells will dominate the stellar surface at one time and produce a temperature inhomogeneity 
as large as 1000 K. Also, some observations were often interpreted in terms of the inhomo- 
geneity. We have already noticed that the different behaviors of water bands (Ryde et al. 
2006) or the measured diameters (Weiner et al. 2003) between the near- and mid-infrared 
can be due to the photospheric inhomogeneity. Also, imaging of nearby supergiants such as 
Betelgeuse (e.g. Young et al. 2000) revealed a presence of bright spots on the hmb-darkened 
disk. So far, however, it is not known yet how these theoretical and observational results can 
be modeled in a unified picture, and this will be an interesting subject to be investigated 
further. 
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5.2. Observations 

The observation relatively free from the difficulty noted in Sect. 5.1, namely an observed 
result generally depends on several model parameters, is the near infrared spectra, which 
suffer little effect of the extended geometry of the outer atmosphere. As we see in Sects. 3 
and 4, the column density can be estimated reasonably well from the Stratoscope data alone, 
and this fact indeed makes the subsequence analysis rather easy. Thus, the Stratoscope data, 
which were observed 40 years ago, are still unique and invaluable as the probes of the outer 
atmosphere. Unfortunately, few observations were done in the near infrared from outside 
the Earth's atmosphere since then. However, the near infrared region can now be observed 
rather well from the dry sites on ground and the near infrared spectra of higher quality can 
be obtained without much difficulty. Such an observations should be quite useful to examine 
the accuracies of the data taken nearly half a century ago. 

Although the spectra in the region long-ward of 2.5 fim. were extensively observed espe- 
cially with the ISO SWS, the H2O 2.7 fim bands are blended with the other molecular lines 
such as of CO and OH while CO and SiO bands include large contributions from the photo- 
sphere. For these reasons, the H2O 1.9 /xm bands are the best probes of the outer molecular 
layers of the early M giant and supergiant stars, and especially the H2O column density can 
be estimated quite well almost independently of other parameters. 

Despite the general difficulty to observe H2O from ground, recent narrow band inter- 
ferometry by Perrin et al. (2005) demonstrated that the observation of the water layer is 
possible with the ground-based interferometer and provided decisive evidence for the extra 
molecular layer in /i Cep. This is a quite encouraging result in that the ground-based in- 
terferometry will already provide fruitful results for the study of the structure of the extra 
molecular layers in cool luminous stars, before the interferometry in space can be realized. 

The H2O ^2 bands in the 6-7 /im region are excellent probes of the geometrical extension 
of the molecular layers, but no longer be useful for estimating the column density because 
they include appreciable emission components from the extended atmospheres and finally 

turn to emission. Extensive observations of this region were first possible with the ISO 
SWS. It is interesting to notice that /x Cep is almost unique in showing the H20z^2 bands in 
emission, except for some Mira variables (Yamamura et al. 1999). Our preliminary survey 
of dozens of red giant stars in the ISO archive revealed no object with H20z/2 bands in 
emission, although the absorption bands are stronger and weaker in the early and late M 
giants, respectively, compared with the predictions based on the photospheric models (Tsuji 
2002b). 

It is to be noted that our modelings of the molsphere are done with the use of the spec- 
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tral and visibility data short-ward of 7.5 //m alone. This has the advantage that the modeling 
can be done almost independently of the effect of the dust components including the newly 
identified alumina dust shell (Verhoelst et al. 2006). Also the absorption components still 
dominate especially in the near infrared region, and this fact makes the diagnosis relatively 
easy. In fact, it is almost impossible to determine the column density, if the emission com- 
ponents dominate as in the longer wavelength region. Thus, the basic parameters of the 
molsphere such as T^^, Ncoi, and i?in are determined from the observations in the shorter 
wavelength region while additional information can be obtained from the observations in the 
longer wavelength region. 

The 10 /xm region can be accessible from ground, and detailed interferometry and spec- 
troscopy realized recently for a Ori provided further constraints as well as new problems on 
the outer atmospheres of red supergiant stars, as discussed in Sect.4.3. It highly desirable 
that similar observation of the spectra and visibility can be extended to the mid-infrared 
region of fi Cep. In the 40 /xm region observed with the ISO SWS, water appears in emission 
in n Cep but only marginally in a Ori (Tsuji 2000b). Our preliminary version of the Model 
A plus the silicate dust shell could reproduce the whole SWS spectrum of fi Cep between 
2.5 and 45 /im reasonably well (Tsuji 2002b). Thus our molsphere model shows overall con- 
sistency with the observed infrared spectrum covering a wide spectral range. However, some 
details such as the relative intensities of some emission lines do not necessarily agree very 
well with the predicted thermal emission spectra, and further detailed analysis including the 
non-LTE effect should be needed. 

The bright red supergiant stars have also been targets of extensive radio observations, 
and one of the highlights may be the VLA mapping of Betelgeuse at 7 mm by Lim et al. 
(1998), who showed that the temperatures over the same height range as the chromosphere 
(e.g. from 2 to ^ 1 R*) are appreciably cooler (correspondingly from 3450 ± 850 K to 
1370 ± 330 K) than the temperatures generally assumed for the chromosphere (e.g. 8000 K). 
Our molsphere may be situated just inside of such radio photosphere^ after a naming by Reid 
& Menten (1997), and covers the region not well resolved with the radio interferometry. A 
detailed modeling of the outer atmosphere of Betelgeuse based on the radio data has been 
done by Harper et al. (2001), who also surveyed a large amount of observations covering 
from UV to cm regions. It seems that inhomogeneity must be considered for all the multi- 
wavelength data to be integrated into a unified picture of the outer atmosphere of red 
supergiant stars. This fact also suggests that the interrelation between the molsphere and 
the stellar chromosphere should be investigated. 
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5.3. Water in Red Supergiant Steirs 

Within the hmitation of the present modehng discussed above, the extra molecular 
layers of a Orionis and Cephei seem to show interesting differences. The molsphere of a 
Orionis is relatively hot (Tex ~ 2250 K) and compact (i?in ~ 1.3 R^) while that of fj, Cephei 
is cooler (Tex ~ 1600 K) and more extended (i?in ~ 2.0 i?*). It is interesting if such a change 
may represent an evolution of the outer atmosphere in the supergiant evolution. 

Probably more advanced stage of the evolution of the outer atmosphere in red super- 
giants may be represented by such objects as VY CMa and S Per. In these objects, infrared 
excess is very large, indicating that the dust envelope should be very thick. The gaseous 
components in the outer atmosphere are more difficult to observe and, even though these su- 
pergiants were extensively observed with ISO (e.g. Tsuji et al. 1998; Harwit et al. 2001), little 
is known yet. On the other hand, these red supergiants are known as the maser sources. Of 
particular interest is that water appears to be quite abundant in the outer space of these red 
supergiants. Recent VLBI observations revealed that water forms masering clouds around 
VY CMa (Imai et al. 1997) as well as S Per (Richards et al. 1999). 

The origin of the masering water clouds may not be known yet. It is an interesting 
possibility that the molsphere in the early M supergiants will develop to the masering water 
clouds in the later M supergiants. But the origin of the molsphere itself is unknown. It may 
also be possible that unknown water supply will provide water to the molsphere as well as 
to the masering clouds. Anyhow some mechanism is needed to enhance the matter density 
in the outer atmosphere around red supergiant stars, either transporting more matter from 
the central star or else accreting some matter from the outer space. Certainly we are still far 
from realizing the recycling of water around evolved stars, and comparative studies of objects 
in the different stages of developing the outer atmosphere will be useful for this purpose. 

6. Concluding Remarks 

The presence of the extra molecular layers in the outer atmosphere of red supergiant 
stars has been anticipated from the infrared spectra for the first time. With the spectra 
alone, however, it was difficult to exclude a possibility that the unexpected spectral features 

may be due to anomalous structures of the stellar photospheres. Recent interferometric 
observations in different spectral regions finally provided decisive evidence for the presence 
of the extra molecular layers outside the stellar photosphere. 

On the other hand, interpretation of the visibility data alone may be by no means unique 
since the present interferometric observations do not yet reconstruct the astronomical image 
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directly. For this reason, simultaneous analysis of the spectral and visibility data should 
be quite essential. Certainly, mult i- wavelength interferometry such as done recently by 
Perrin et al. (2005) already includes some spectroscopic information in itself, and future 
spatial interferometry with a higher spectral resolution will hopefully involve all the necessary 
spectroscopic information. 

Thanks to the fine achievements, both in the spectroscopy and interferometry, we are 
convinced of the presence of the extra molecular layers outside the stellar photosphere. Now, 
an important problem is how to understand the presence of such a rather warm and dense 
molecular layer outside the photosphere. This problem may be related to that of the stellar 
chromosphere, for which detailed modelings have been done but its origin may be by no 
means clear yet. Probably, the problem of the origin of the molsphere may be as difficult 
as that of the chromosphere, and will require comprehensive understanding of whole the 
outer atmosphere. We hope that more attentions, both observational and theoretical, could 
be directed to this problem of the extra molecular layer, or the molsphere, in view of the 
convincing confirmation for its existence in the outer atmospheres of red supergiant stars. 

I would like to thank an anonymous referee for careful reading of the text and for 
helpful suggestions resulting in improving the text. I also thank Dr. N. Rydc for making 
available the preprint of his recent paper, which stimulated the present study, at least partly. 
My thanks are also due to Dr. K. Ohnaka for reading the draft of this paper with useful 
comments. Data analysis has been done in part on the general common-use computer system 
at the Astronomical Data Analysis Center (ADAC) of NAOJ. This work is supported by the 
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Table 1: Basic Stellar Parameters 



item 


assumed value 




15.0 




650 




3600 K (a Ori) , 3800 K {^l Cep) 


Cmicro 


5.0 kms~^ 


Cmacro 


10.0 kms"^ 


log Ac /An 


-4.0 


log AnMh 


-3.5 


log Ao/A^ 


-3.5 



Table 2: Empirical Models of the Molsphere around /j, Cephei 



Models Model A Model B Model C Model D'^ Model D* 



Teff(K) 


3600 


3800 


3800 


3800 


3800 


/ Rq 


650 


650 


650 


650 


650 


Tex(K) 


1500 


1500 


1600 


2700 


2700 


Rin/ Rq 


1306 


1438 


1310 


851 


851 


Rout/ Rq 


1699 


1700 


1723 


858 


858 


iVeol(H2O)/1020 


3.0 


3.0 


3.0 


0.5 - 33.0 


2.8 


iVeol(CO)/1020 


3.0 


3.0 


3.0 


33.0 


93.0 



"Based on the parameters suggested by Perrin et al.(2005) 
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Table 3: Optical Thickness and Column Density 



wavelength (/im) 




molecule 


K(cm^)^ 




A^eoi (cm-2) 


2.03 


0.22 


H2O 


1.6 X 10- 


-21 


1.4 X 1020 


2.15 


0.02 


H2O 


4.5 X 10" 


-22 


4.4 X 10^9 


2.22 


0.07 


H2O 


2.5 X 10" 


-22 


2.8 X 10^0 


2.39 


2.62^ 


H2O 


8.0 X 10" 


-22 


3.3 X 10^1 




i.sr 


CO 


4.0 X 10- 


-22 


3.3 X 10^^ 



"Perrin et al.(2005) 
''read from Fig.5 

''The total optical depth of r = 3.93 by Pcrrin ct al.(2005) is divided into t(H20) = 2.62 and r(CO) = 1.31 
so that the column densities of H2O and CO are the same. 



Table 4: Empirical Models of the Molsphere around a Orionis 



Models Model E Model F Model G Model H 



Teff(K) 


3600 


3600 


3600 


3600 




650 


650 


650 


650 


Tex(K) 


1500 


1750 


2000 


2250 


Rin/ Rq 


1107 


980 


914 


850 


-Rout/ R(3 


1706 


1766 


1803 


1849 


iVeoi(H2O)/1020cm-2 


1.0 


1.0 


1.0 


1.0 


7V,„i(CO)/1020cm-2 


1.0 


1.0 


1.0 


1.0 
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4000 
T (K) 

3000 
2000 
1000 



-10 -5 ' 

log Tq 

Fig. 1. — The temperature T, radius R, and logarithms of the specific density logp of the 
LTE classical model photosphere in radiative and hydrostatic equilibrium (M = 15 Mq, 

i?* — 650 i?0, Tcfr = 3600 K ) plotted against logro (tq is the optical depth in the continuum 
opacity at 0.81 /xm). The solid lines are the case that the integration started at logro = -6-0 
and this is what usually done for photospheric models. The dashed lines show the case that 
the integration started at logro = -14.0. The model is still within the generalized Eddington 
limit and remains stable out to ~ 1700 i?© ~ 2.6 i?* at least formally. 






A(/xm) 

Fig. 2. — (a) The predicted band averaged visibilities (solid lines) based on the Model A for 
molsphere + photosphere {T^s = 3600 K, Tex = 1500 K, ^ 2.0 R^, A^'coKHzO) = A/'coi(CO) 
= 3.0 X 10^° cm~^) are compared with the observed results for fi Cep (filled symbols) in the 4 
narrow bands by Perrin et al. (2005). (b) The predicted near infrared spectrum based on the 
Model A is compared with the spectrum of /i Cep observed with the Stratoscope (resolution 
R ^ 200) and corrected for the interstellar reddening with A^j = 1.5 mag. The predicted 
spectrum of the photosphere, which shows only CN and CO bands but not H2O, is also 
shown, (c) The predicted spectrum based on the Model A is compared with the spectrum of 
fi Cep observed with the ISO SWS (resolution R ^ 1600) and corrected for the interstellar 
reddening with = 1.5 mag. The predicted spectrum of the photosphere appears at the 
lower left corner. 






A(//in) 

Fig. 3.— (a) The same as Fig.2a, but for the Model B {T^s = 3800 K, = 1500 K, 
i?in ^ 2.2 R^, NcoiiW) = A/'coi(CO) = 3.0 x lO^Ocm-^). (b) The same as Fig.2b, but for 
the Model B. (c) The same as Fig.2c, but for the Model B. 
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A(/zm) 

Fig. 4.— (a) The same as Fig.2a, but for the Model C (Teg = 3800 K T^^ = 1600 K, 
^ 2.0 N^oiiW) =N^oi{CO) = 3.0 X 1020cm-2). (b) The same as Fig.2b, but for the 
Model C. (c) The same as Fig.2c, but for the Model C. 
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h.6 1.8 2 2.2 2.4 

A(/U.m) 



Fig. 5. — a) Absorption cross-section (cm^ molecule"^) of H2O at T = 2700K. The black 
line is by the high resolution {R ^ 50000) and the white line is the straight mean opacity 
smeared out with the resolution it! ^ 1600. b) The same as for a) but for CO. 
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Spatial Frequency (cycles/arcsec) 



Fig. 6. — The predicted band averaged visibilities based on the monochromatic visibihties 
are shown by the sohd hnes for the Model D (Teg = 3800 K, Tex = 2700 K, R^^ ^ 1.3 R^, and 
the optical thickness for the 4 narrow band regions given in Table 3). It is confirmed that 
the results for the K215, K222, and i^239 bands based on the column densities estimated 
with the cross-sections of the contributing molecular bands (Fig. 5) reproduce the predicted 
visibilities based on the mean optical thicknesses (Fig.l of Perrin et al.(2005)). The results 
also agree with the observed visibilities for n Cep by Perrin et al (2005) shown by the filled 
symbols. But the result for the K203 band for logA'coi = 1-4 x 10^°cm~^ in Table 3 (solid 
line labelled with 1) does not reproduce the result of Perrin et al.(2005). Instead we find that 
reasonable agreements with the predicted as well as observed values by Perrin et al.(2005) 
can be obtained for logAi'coi = 4.2 x 10^° cm~^ (dotted line). 
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A(/i,m) 

Fig. 7.— (a) The same as Fig.2a, but for the Model D* (Teg = 3800 K, Tex = 2700 K, 
^ 1.3/2,, N^oiiW) = 2.8 X 1020cm-2, Ar,oi(CO) = 9.3 x lO^^cm-^). (b) The same as 
Fig.2b, but for the Model D*. (c) The same as Fig.2c, but for the Model D*. 
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A(/im) 

Fig. 8. — (a) Predicted spectra of water in the 12 fim region (spectral resolution of 0.005 cm~^ 
and broadened with the macro-turbulence of lOkms"^) based on the LTE classical model 
photospheres with the parameters in Table 1 except for Tefj's, which are 1: 3600 K, 2: 3400 K, 
3: 3200 K, 4: 3000 K, and 5: 2800 K. Some lines of water can be seen already at T^s ~ 3600 K 
and stronger in the cooler models, (b) The same as (a), but for the near infrared region 
(spectral resolution of 0.1 cm~^ and convolved with the slit function of FWHM = 1500 
kms~^). Since the /-values of H2O lines are about 1-2 orders of magnitude smaller in this 
region compared with those in the the 12 /im region, water bands can be seen only in the 
models with T^s < 3000 K. 
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Fig. 9. — (a) The predicted near infrared spectrum based on the Model E (Tgx = 1500 K, 
i?in ~ 1.7 i?^,, Ncoi = 1-0 X 10^°cm~^) is compared with the spectrum of a Ori observed 
with Stratoscope (resolution R ^ 200) and corrected for the interstellar reddening with 
= 0.5 mag. The predicted spectrum of the photosphere is also shown, (b) The predicted 
spectrum of the H2O 1^2 band region based on the Model E is compared with the spectrum 
of a Ori observed with the ISO SWS (resolution R ^ 1600) and corrected for the interstellar 
reddening with = 0.5 mag. The predicted spectrum of the photosphere appears below 
that for molsphere + photosphere. 
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Fig. 10.— (a) The same as Fig.9a, but for the Model F (Tg^ = 1750 K, R^^ ^ 1.5 R^, N^^i 
= 1.0 X 10^°cm~^). The dotted hne is for the case of A^'coi = 2.0 x 10^^ cm^^ suggested by 
Verhoelst (2006). (b) The same as Fig.9b, but for the Model F. 
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Fig. 11.— (a) The same as Fig.9a, but for the Model G (Tex = 2000 K, ^ 1.4 A^'^oi 
= 1.0 X 10^° cm"^). The dotted hne is for the case of A^coi = 2.0 x 10^° cm"^ suggested by 
Ohnaka (2004). (b) The same as Fig.9b, but for the Model G. 
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Fig. 12.— (a) The same as Fig.9a, but for the Model H (Tex = 2250 K, ^ 1.3 R^, N^^i = 
1.0 X 10^0 (b) The same as Fig.9b, but for the Model H. 
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Fig. 13. — The predicted visibilities for the K band based on the molsphere + photosphere 
models (Models E- H) (solid lines) are compared with the observed ones in the K bands (filled 
circles) by Perrin et al. (2004). Also, the predicted band averaged visibility curve (solid line) 
based on the classical photospheric model (M = ISM©, = 650 Rq, Tes = 3600 K, and 
other parameters in Table 1) is shown for comparison. 



-47- 




Fig. 14. — The predicted spectrum around SiO fundamental bands based on the Model H 
is compared with the spectrum of a Ori observed with the ISO SWS (resolution R ^ 1600) 
and corrected for the interstellar reddening with Ay = 0.5 mag. The predicted spectrum of 
the photosphere appears below that for the Model H. 



